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ABSTRACT 

We have carried out a spatial-kinematic study of three proto-planetary nebu- 
lae, IRAS 16594—4656, Hen 3-401, and Rob 22. High-resolution H 2 images were 
obtained with NICMOS on the HST and high-resolution spectra were obtained 
with the Phoenix spectrograph on Gemini-South. IRAS 16594—4656 shows a 
"peanut- shaped" bipolar structure with H 2 emission from the walls and from 
two pairs of more distant, point-symmetric faint blobs. The velocity structure 
shows the polar axis to be in the plane of the sky, contrary to the impression given 
by the more complex visual image and the visibility of the central star, with an 
ellipsoidal velocity structure. Hen 3-401 shows the H 2 emission coming from the 
walls of the very elongated, open-ended lobes seen in visible light, along with a 
possible small disk around the star. The bipolar lobes appear to be tilted 10—15° 
with respect to the plane of the sky and their kinematics display a Hubble-like 
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flow. In Rob 22, the H2 appears in the form of an "S" -shape, approximately 
tracing out the similar pattern seen in the visible. H 2 is especially seen at the 
ends of the lobes and at two opposite regions close to the unseen central star. 
The axis of the lobes is nearly in the plane of the sky. Expansion ages of the 
lobes are calculated to be ~1600 yr (IRAS 16594-4656), ~1100 yr (Hen 3-401), 
and ~640 yr (Rob 22), based upon approximate distances. 

Subject headings: circumstellar matter — infrared: stars — ISM: individual (Hen 
3-401, IRAS 16594—4656, Rob 22) — planetary nebulae: general — stars: mass 
loss — stars: post-AGB 



INTRODUCTION 



Proto-planetary nebulae (PPNs) are objects in transition between the asymptotic giant 
branch (AGB) and planetary nebula (PN) stages of stellar evolution. Studies over the past 
decade, particularly imaging studies with the Hubble Space Telesco pe (HST), have show n 
that this transition stage is key to unde rstanding the shaping of PNs (IBalick fc Frankl|2002l ). 
In the interacting stellar winds model (Kwok Ill982l ). which was later generalized to include 
an equatorial density gradient (jBalicki 1 198 71 ). it was assumed that a fast wind from the 
central star interacted with the detached, slowly-moving remnant AGB envelope to shape 
the nebula. HST i mages of PPNs show that many and perhaps most of them display a basi c 
bipolar structure JUeta et al.l I2OO0I : ISu et al.l lioOll : ISahai et alil2007l : Ishdmiak et "all l2008h . 
These are shapes that appear to be further dev eloped in the PN stage. In addition, a point 
symmetry is often seen; ISahai fc Traugerl (119981 ) attribute this to collimated jets, which may 
be episodic and change their direction, and which they advocate as the main shaping agents 
of PNs. 



Studying H 2 in PPNs is one of the best ways to probe the presence of a fast, per- 
haps collimated, post-AGB wind and its effect in shaping the nebula. H 2 is the main con- 
stituent of the detached circumstellar envelope around a PPN. While a fast wind (V>100 
km s _1 ) has the energy to dissociate the H 2 on direct impact, it can also produce shocks 
in the medium that move at a slower speed and collisionally exci te the H 2 . Recent su r- 



veys have shown the presence of shock-excited H 2 in bipolar PPNs (iGarcfa-Hernandez et al. 



20021 : iKelly &: Hrivnakl 120051 ). Det ailed high-resolution H 2 studies have been published of 



half dozen PPNs ( 


Sahai et al. 


19981: 


Kast 


ner et al. 


2001 


Hrivnak et al. 


20061: 


Van de Steene et al. 


2008|). 



2003 



In this paper we present the results of a detailed H 2 emission study of three bipo- 



- 3- 



lar PPNs, IRAS 16594-4656 ("Water Lily Nebula"), Hen 3-401 (IRAS 10178-5958), 
and Rob 22 (IRAS 10197-5750). All three have relatively hot central stars for PPNs, 
B and A spectral types, and all possess large infrared ex cesses due to circumstellar dust 
( IParthasarathy &: Pottaschlll989l ; IGarcia-Lario et al.lll999a( ). All three are or would be clas- 
sified as DUPL EX ("DUst Prom inent, L ongitudinally- Extende d" ) nebulae in the classifica- 



tion system of lUeta et al.l (120001 . see also IShdmiak et al.l (120081 )). 



Rob 22 has a nebula that looks like a pair of butterfly wings. At high spatial resolution 
the lobes show a filamentary structure in visible light. In addition to bipola r structure, 
some p oint symmetry is apparent in the lobes. A large halo extends out to 25" ( ISahai et al. 
1999al ). The nebula appears to be viewed essentially edge on, with a dark lane o bscuring the 
star in visible light. In the more detailed morphological classification system of ISahai et al. 
( 120071 ) . this object is classified as Bcw,ml,ps(s),h(e): bipolar nebula with closed lobes, cen- 
tral obscuring waist, minor lobes present, point-symmetric shape, and enlongat ed halo. Th e 
spectral type is A2 Ie, based on light reflected off the lobes of the nebula (j Allen! 1 19781 ). 
The circumstellar envelope appears to have a d ual chemistry, with an oxygen-r ich region 
evidenced by OH maser e mission (OH284. 18- 0.79; lAUen. Hyland. &: Caswelllll980l ) and crys- 
talline sil icate emission (IMolster et al.ll2002l ) and a carbon-rich region evidenced by PAH 



emission (IMolster et al.lll997l ) 



The nebula of Hen 3-401 possesses long, narrow lobes. It seems to be viewed at a 
slightly larger angle with res pect to the plane o f the sky, with the central star appear- 



ing faintly between the lobes (ISahai et al. 



1999b! ) . The morphological classification of this 



object is Bow*(0.6),sk: bipolar nebula with open lobes, central ob scuring waist with star 
evident (at 0.6 //m), with a skirt-like structure (ISahai et al.l 120071 ). The spectral type of 
the star is Be, and the spectrum shows stron g Balmer emission lines and lots of per mitted 
and forbidden low-excitation emission lines (lAllenl 1 19781 ; IGarcia-Lario et al.l Il999bt). The 



circumstellar envelope appears to be carbon- rich, since it shows CO emission (ILoup et al. 



1990 



Buiarraba 



OH (ISilva et al. 



&; Bachillerlll99ll ) and PAH emission (IParthasarathy et a 



19931 ) emission or crystalline silicates (IParthasarathy et al 



200lh . 



200l[) but not 



In IRAS 16594—4656, the star is clearly seen and relatively bright compared to the 
nebula, giving the impression that the nebula is at a larger orientation angle. The nebula 
has the appearance of a basic bipolar structure with a pronounced point symmetry, consist- 
ing of three pairs of oppositely-directed, slightly curved, thin lobes. Also seen is a s maller 
elliptical structure oriented from northwest to southeast (IHrivnak. Kwok. fc Sulll999l ). The 
morphological classification of this object is Mcw*(0.6),an,ps(m,an),h(a): multipolar neb- 
ula with closed lobes, central obscuring waist with star evident (at 0.6 /mi), ansae present, 
point symmetry with two or more pairs of diametrically opposed lobes and with ansae, and 
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a halo wit h centrosymmetric arcs (Sahai et al.l 120071 ). Its spectrum shows Balmer emis- 
sion lines ( Garcfa-Lario et al. 1999af) and has been classifi ed as B7 based on the Ha line 



( IVan de Steene. Wood. &: van Hoofll2000l ). iReyniersI (120021 ). on the basis of high- resolution 
spectra, fit the spectrum with T eff =14,000 K and log<7=2.1, slightly hotter than (5 Ori (B8 I); 
this indicates a spectral type of B5-7 Ie. He also observed many emission lines. The neb- 
ula appears to be carbon rich, showing PAH features and displaying the 21 jum and 3 
ium features seen in carbon -rich evolved stars (jGarcia-Lario et al.l Il999al ; IVolk et al.ll2002l ). 
Van de Steene et al.l (120081 ) recently published a detailed kinematic and morphologic study 
of this object using H2 emission, and we will refer to their results in the present study. 



We begin with an examination of the high-resolut ion H-? images ( see ISahai et al.l (120001 ) 
for preliminary images of Rob 22 and Hen 3—401 and iHrivnak et al.l (120041 ) for preliminary 
images of IRAS 16594—4656), then present new high-resolution, long-slit H2 spectra, and 
then follow with a discussion of the kinematics. From these, we are able to determine 
properties of the collimated winds in these three PPNs and how they have shaped the 
surrounding nebulae. 



2. HIGH-RESOLUTION H 2 IMAGING OBSERVATIONS 

High-resolution near-infrared imaging observations of the three PPNs were obtained 
with the Hubble Space Telescope (HST) and the Near-Infrared Camera and Multiobject Spec- 
trometer (NICMOS) camera 2 (NIC2). The camera has a resolution of 0'.'076 x 0'.'075 pixel -1 
and a field of view of 19'.' 5 x 19'.'3. Observations were made with a narrow-band H 2 filter 
(F212N) and also with an adjacent continuum filter (F215N) for the purpose of distinguishing 
the H 2 1—0 S(l) emission from the scattered starlight. 

Hen 3-401 and Rob 22 were observed on 1998 under General Observer program ID No. 
7840 (PI: Kwok) using the original HST NICMOS camera. Observations were also made 
with the F160W (H) and F222M (K) filters for both of these objects. IRAS 16594-4656 
was observed in 2002 under General Observer program ID No. 9366 (PI: Hrivnak) using the 
refurbished HST NICMOS camera. Observations were also made with the F110W (J) filter. 
The observing log is listed in Table HJ 

The data were taken in the multiple non-destructive mode (MULTIACCUM) to obtain 
a high dynamic range without saturating the detector. Three sets of observations were made 
through each filter at each position on the array, and each target was imaged at three (Hen 
3-401, Rob 22) or four (IRAS 16594-4656) different positions of the NIC2 array using a 
predefined dither pattern. 
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The data were processed using the standard NICMOS calnica pipeline, which applied 
bias subtraction, dark subtraction, flat-fielding, and cosmic ray removal. The more recent 
data for IRAS 16594—4656 were re-calibrated using the HST archival On-the-Fly Repro- 
cessing system, which provides recalibration using the most up-to-date reference files. The 
individual files were then corrected for the "pedestal effect" using IRAF/STSDAS pedsub 
task to flattening the sky value in the four quadrants of the NIC2 array. We then com- 
bined the twelve images with each filter using the I RAF STSDAS/Dither pa ckage with the 



parameters suggested in the HST Dither Handbook ( jKoekemoer et al.ll2002l ). The data for 



Hen3— 401 and Rob 22, which were taken and reduced ea rlier, were reduc ed with the same 



philosophy but in a more manual fashion, as described by ISu et al.l (120031 ) 



The reduced F212N and F215N images were then flux calibrated using the appropriate 
calibration values for each filter before and after the NICMOS refurbishing. These calibration 
values are listed in Table [2j The individual F212N and F215N images for the three PPNs, 
with the thermal backgrounds removed, are displayed in Figure [TJ 

For IRAS 16594—4656, the relatively-bright central star is surrounded by a distinct 
peanut-shaped nebula in the F212N image, while only some smaller and fainter nebulosity 
surrounds the star in the F215N image. Thus one can immediately see that most of the 
nebular light in the F212N image is due to the H2 emission rather than scattered starlight. 
One can also see faint emission from two of the pairs of longer, thin lobes, with brighter 
emission ("ansae") at their ends. For Rob 22, the F212N and F215N images appear rather 
similar. They each show, rather than a central star, emission coming from two closely-spaced 
regions separated by a dark equatorial band. These bright regions may represent light from 
the central star scattered toward us from above and below an obscuring disk or torus. The 
two bright regions lie roughly along the line that connects the approximately point-symmetric 
outer ends of the two lobes. The dark space between the two bright regions might thus be 
due to the torus. The nebula appears to be a little brighter in the F212N image than in the 
F215N image, especially in the outer regions, indicating that in the F212N image the nebula 
is seen mostly in scattered light but with some contribution from H2 emission. For Hen 
3—401, the nebula is seen more clearly and is larger in the F212N image than in the F215 
image, indicating that in the F212N image, it is seen in both scattered light and H 2 emission. 
The star is bright in each image. 



3. H 2 IMAGES 



To produce images of these three PPNs showing only H2 emission requires the removal 
of the scattered light from the F212N image, which can be accomplished through the sub- 
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traction of an appropriately scaled image of the adjacent continuum (F215N image). The 
scaling is necessary to correct for differences in the source continuum in the two bandpasses. 
Differences due to the filter profiles and detector respon se were corrected by the flux cali- 



bration. In a similar previous study (IHrivnak et al.ll2006l ). we examined the question of how 



best to scale the F215N continuum image for subtraction from the F212N image to remove 
the scattered continuum light. The object of that study, IRAS 17150—3224, had no other 
emission lines in these two bandpasses, and thus a simple correction for the continuum levels 
was appropriate. In these cases, we also need to be concerned about other emission features 
which may be contributing to the flux in the two filters. 

In the following sub-sections we discuss the details of the continuum removal and the 
resulting H 2 image for each object. 



3.1. IRAS 16594-4656 



A spectrum for IRAS 16594—465 6 in the 2 /xm region (2.10—2.27 /mi) has been published 
by lVan de Steene Sz van Hood (b()03h . It was obtained with a long slit l'.'O wide and oriented 



E-W. The spectrum shows a relatively strong H2 feature at 2.12 /im with a much weaker Br7 
emission line. Since the Br7 is at a wavelength that places it at the edge of a wing of the 
F215N filter profile, its contribution is insignificant in this image. This is illustrated in Figure 
[21 which shows the spectrum of the object plotted together with the filter profiles. Thus, 
for the continuum removal, we convolved a linear fit to the observed spectral continuum 
(with the H 2 emission removed) with the filter profiles and the detector response using the 
STSDAS/SYNPHOT calcphot task. This resulted in a scale factor of 1.016 that we used to 
scale the F215N image prior to subtraction of the flux-calibrated images. 

The resultant H2 image of IRAS 16594—4656 is displayed in Figure [3^. It shows a limb- 
brightened, double-lobed bipolar structure with the west lobe larger and brighter (ratio of 
1.6) than the east lobe. Approximate sizes (lOer level) are 2'.'6 by 2'.'8 for the west lobe and 
2'.'2 by 2'.'5 for the east lobe. In addition, one can see two pairs of more distant H 2 -emitting 
ansae, each pair collinear on a line passing through the star. The ansae to the SW are 
brighter and farther from the star than those to the NE, with distances of 6'.'2 and 5'.'4 for 
the ansae at PA=33° and of 5'.'8 and 4'.'7 for the those at PA=54°. These ansae are found 



at the ends of the thin lobes seen in the visible light images (IHrivnak. Kwok. fc Sulll999 



Su et al.l l200ll ). and will be discussed together with the visible image in Section [5j The 
emission at the position of the star appears to approximately cancel out, indicating that 
there is little, if any, H 2 emission from the star. 
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3.2. Hen 3-401 

Spectra of the 2.12 ixm H 2 line and the Br7 line of Hen 3-401 have been published by 



Garcfa-Hernandez et al.l (120021 ) . They observed this object twice with a 4'.'5-wide slit oriented 
at P. A. = 70° and 90°, approximately along the bipolar axis, and they found the H 2 emission 
to be clearly extended. The H 2 line is stronger than the Br7 line by about a factor of 1.5, 
and since the latter lies in the wing of the F215N filter, its contribution to the flux in the 
F215N image is very small. This is illustrated in Figure [2j Neglecting the contribution of 
the line results in a scale factor of 0.994 (the spectrum is essentially flat in this region). 

The resulting H 2 image of Hen 3-401 is shown in Figure [3b. H 2 emission extends along 
the two almost cylindrical, slightly barrel-shaped lobes, each with a size of 2'.'5 x 10". The 
east lobe is brighter, with a ratio of 1.2. The two lobes are limb-brightened and appear to 
be brightest 2V 5— 5'.' 3 (E) and 1'.'5— 4'.' 3 (W) from the star. The ends of the lobes are open. 
Little, if any, of the H 2 emission appears to come from the central star[j] 



3.3. Rob 22 



Rob 22 was also observed spectroscopically in this region by iGarcfa-Hernandez et al. 



( 120021 ). They observed with the slit oriented at P.A. = 90°, in this case almost perpendicular 
to the bipolar lobes. With a slit width of 4'.'5, they likely included most of the nebula but 
missed some of the H 2 emission from the very ends of the lobes. Their spectra show the Br7 
line to be stronger than the H 2 2.12 fim H 2 line by about a factor of 2.5, but since it still 
falls on the extreme edge of the filter profile, its contribution will be small. Nevertheless, we 
examined the scaling of the F215N image with and without consideration of this emission 
line. Neglecting the contribution of the line results in a scale factor of 0.979, while including 
it results in a scale factor of 0.967. We will use the scale factor with the line included, but 
note that the difference is only at the 1% level. Note, however, that the difference can be 
larger in localized regions, since the Br7 and H 2 emissions likely have different distributions. 

The resulting H 2 image of Rob 22 is shown in Figure [3fc. The nebula appears to have 
somewhat of a point-symmetric "S" shape, with a bright clump in the inner part of each 
lobe near the position of the unseen central star (clump to the NE is brighter) and a ridge of 
H 2 along the the end of each lobe (especially the S lobe). The north lobe is brighter, with a 
ratio of 1.2. Much of the extended structure within a few arcsec of the star disappears after 



1 Garcia-Lario et al.l (|1999bT ) previously displayed the HST H2 data, using a scale factor based upon the 



field stars, and obtained a similarly appearing H2 image. 



- 8- 



continuum subtraction, indicating that the features at the bases of the polar lobes nearest 
the star are due to scattered light from a dusty flared disk or torus, not H 2 emission. 



4. HIGH-RESOLUTION H 2 SPECTRA 

4.1. Observations and Reduction 

High dispersion, long-slit spectral observations were made of the three sources in the 
H 2 v=l— S(l) lin e at 2.1218 /mi using the near-infrared echelle spectrograph Phoenix 



( IHinkle et al.lll998l ) on the 8-m Gemini-South telescope on Cerro Pachon, Chile. Phoenix 
has a 1024x256 InSb detector with a pixel scale of 0'.'085 x 1.26 km s _1 at a wavelength 
of ~2 /mi. The slit width was (X'34, providing an effective velocity resolution of 6 km s _1 
(R~50,000). The spectrum covers a range of ~1300 km s^ 1 (0.0086 /mi) centered on this 
H 2 line. The area of the sky projected by the long slit onto the detector is ~13'.'4 in length, so 
for the observations of IRAS 16594—4656 and Hen 3-401, more than one exposure was taken 
moving along the slit, yielding a longer effective slit to cover each target. The observations 
were made in 2003 under three queue observing programs (GS-2002B-Q-53, GS-2003A-Q-27, 
GS-2003B-Q-41). Image quality on the different nights varied from 0'.'7 to l'.'O as measured 
in the R band. 

Several spectra were obtained for each of the PPNs at various positions. For each 
object, one spectrum was taken with the slit oriented through the star along the major axis 
of the lobes. Additional spectra were taken at various offsets and orientations. For IRAS 
16594—4656, two of the slit positions were oriented to pass though the star and the pairs of 
H 2 clumps. An observing log listing the slit positions is given in Table [3j and the orientations 
of the slit for these various observations are shown in Figure [31 

The weather for these observations was mostly clear, although some patchy clouds were 
present during observations of Rob 22 and IRAS 16594—4656. Sky subtraction was ac- 
complished by obtaining one exposure of blank sky near each source immediately before or 
after a pair of exposures of the target. Spectral flat fields were obtained using an internal 
quartz emission lamp. For each target, a bright standard star was observed on the same 
night with the same grating setting in order to correct for telluric absorption and for flux 
calibration. Numerous telluric lines were also used for wavelength calibration, adopting the 
telluric spectrum available from NOAO0. We calculated the apparent velocities adopting a 
rest wavelength of 21218.356 A for the H 2 v=l-0 S(l) line (Bragg, Brault, & Smith 1982), 



2 ftp : / /ftp . noao.edu/catalogs / arcturusatlas /ir / 
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and these velocities were then corrected to a heliocentric reference frame. Uncertainty in the 
resulting velocities is roughly ±1 km s _1 , dominated by scatter in the dispersion solution for 
the numerous telluric lines across the observed wavelength range. 

4.2. Position- Velocity Diagrams 

4.2.1. IRAS 16594-4656 

Spectra of IRAS 16594—4656 were obtained at five different positions on the nebula, as 
illustrated in Figure 3. The resulting position-velocity (PV) diagrams are shown in Figure 
HI In all cases the brightest emission displays a closed ellipse, but with different sizes and 
strong variations in intensity around the ellipse. In slit position a (PA=72°), which is along 
the major bipolar axis, and slit position b (PA=52°), which is rotated 20° from the major 
axis, the inner ellipses seem to be tilted slightly, with the W-SW lobe slightly blue-shifted 
compared to the E-NE lobe. However, this is a small effect, and the polar axis appears to 
lie very close to the plane of the sky. 

Slit positions c (PA=33°) and d (PA=345°) also pass through the star, but are oriented 
closer to the minor axis (slit position d is along the minor axis), and the size of the bright 
inner ellipse is correspondingly smaller. Slit position e (PA=345°) is parallel to the minor 
axis but offset from the star through the larger W lobe. The PV diagram at this position is 
consistent with a slice through a shell expanding with a velocity of ~8-9 km s _1 and with 
the N side brighter. 

Intensity contours are shown in Figure H] to display the faintest emission levels. Slit 
position a shows that the H 2 emission extends ±5" along the major axis. Slit positions b 
and c pass through the clumps at the ends of extended arms in the nebula, and in both 
spectra these clumps are seen in H 2 emission. In slit position b, the centroid velocities of 
the two clumps differ by 10 km s^ 1 (—36 to —26), with each about 5 km s _1 from the mean 
velocity of ~— 30 km s" 1 . In slit position c, the velocity of the two clumps is about the same. 
These velocities can be explained by a model in which the clumps are ejected by rotating 
jets. The clumps at slit position b appear closer to the star than those at slit position c; this 
can be a result of a projection effect and/or temporal effect, assuming that they are moving 
outward from the central star at the same speeds. If the clumps in each point-symmetric 
pair began moving from the location of the star at the same time, then they should be at 
the same distance from the star. However, in both cases the brighter SW clumps appear 
farther from the star than the NE clumps, so the assumption of constant velocity for the 
clumps may not be correct. Alternatively, the clumps may simply be non- uniform brightness 
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enhancements at the ends of the filaments. 

The bright ellipses are consistent with limb-brightened lobes. When allowance is made 
for the resolution of the spectra, it appears that the interiors of the lobes are devoid of 
H2 emission. In the four slit positions that include the star (Fig. IH panels a — d), the blue- 
shifted side of the inner ellipse is brighter. This may be due to extinction within the lobes 
or intrinsically brighter H 2 emission from the side facing us. The elliptical morphologies of 
the PV diagrams for the four slit positions through the star, especially the one along the 
major axis, are consistent with an expanding ellipsoid. Thus they differ from the "peanut- 
shape" PV diagram that one might expect from two bipolar expanding bubbles. Thus the 
expansion along the minor axis does not appear to be significantly retarded by a constraining 
higher- density torus. 

We investigat ed the velocity field in more detail using the morpho-kinematic modeling 



program SHAPE (ISteffen fc L6peal2006l ). In this program, the morphology of the nebula 
and the velocity field can be specified and a resulting PV diagram produced. The shape 
of the nebula was assumed to be a three-dimensional symmetrical dumbbell based on the 
H 2 (cross-sectional) image. It was then assumed that the H 2 emissivity uniformly followed 
the contours of the nebula. We began by investigating a velocity field in which the velocity 
increased linearly with distance (so-called "Hubble flow"), which is seen in some PNs and 
is seen in Hen 3—401 (Section 4.2). This produced a PV diagram that had a peanut shape, 
contrary to the elliptical shape seen (see Figure left and right panels). When we assumed 
on the other hand a constant velocity (9 km s _1 ), the PV diagram agreed more closely with 
the observed ellipse but was slightly diamond-shaped. A very good fit was obtained by the 
addition of a slight Hubble flow, with V = (8 + 2.4r/l") km s _1 , as shown in the middle 
panels of Figure [5j We determine, based on adding an inclination to the model, that the tilt 
of the polar axis is no more than ~10° from the plane of the sky. 

Slices in velocity across the P-V diagrams at position along the slit = 0" (along the line 
of sight through the star) with width C'77 are displayed in Figure [6j These show a uniform 
expansion of the shell with an expansion velocity of V exp = 9 km s _1 . The velocity center is 
at —30 km s _1 . The expansion velocity of the lobe at 0'.'9 west of the star is similar; this is 
consistent with uniform expansion of an elliptical bubble. 

These observations reveal a systemic heliocentric radial velocity of Vh e i(H 2 ) = —30 km 
s _1 , which translates to Vlsr = ~ 26 km s _1 , with V cxp = 9 km s _1 . These velocities can be 
compared with the mill imeter- wave CO J—l— emission line of Vlsr(CO) = —26 km s _1 



with V oxp = 16 km s _1 ( jLoup et al.lll990l ). The CO expansion represents the velocity of the 
outer remnant of the AGB mass loss expanding into the ambient ISM, while the H 2 expansion 
represents the velocity of the shock front at the edge of the inner bipolar cavity along the 
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line-of-site direction, which is essentially perpendicular to the polar axis. The difference 
between these two expansion velocities appears to reflect the two different kinematic regions. 
Heliocentric radial velocity measurements of individual lines in the visible spectrum showed 
a range of —10 to —40 km s -1 , with aver age values of ^ —20 km s _1 for the absorption lines 
and ~— 32 km s _1 for the emission lines (IReyniersI 120021 ) . 



4.2.2. Hen 3-401 

Spectra of Hen 3-401 were obtained at four different positions on the nebula. The 
resulting position-velocity diagrams for the three oriented along the bipolar axis are shown 
in Figure [7J (The fourth one, passing through the star but oriented perpendicular to the 
axis, showed no extended structure, only H 2 emission at the position of the star with the 
systemic velocity). In each of these three slit positions, we see a pair of convex arcs from 
each lobe. They are not closed at the ends, indicating open lobes out to at least ±7" from 
the star. At the two offset positions, some emission appears to partially fill the interior of 
the lobes; this may simply be contamination from the very edges of the lobes due to seeing. 

The bipolar axis is obviously tilted from the plane of the sky, with the W lobe moving 
toward us, and with velocities along the lobes that are proportional to distance from the 
star (Hubble flow). The pairs of convex arcs represent emission from the thin front and back 
walls of slightly barrel-shaped lobes. The difference between the pair in each lobe is about 
30 km s _1 for the slit passing through the star (Fig. [7b) but about 23 km s _1 in the two 
offset positions; this is consistent with seeing the radial expansion along the bipolar axis 
but only a projected component of it in the offset positions. If the expansion really is a 
Hubble flow, then the relationship between observed expansion velocities of the front and 
back sides of the lobes in Figure [7b compared to their apparent spatial size can constrain the 
dynamical age of the nebula for a given distance. If we assume axial symmetry — so that 
the lobes have roughly the same length-to- width ratio as in H 2 images (this seems justifiable 
based in the similar morphology in Figures [7b and [T2"j) — then one can stretch the PV plot 
in Figure [7b to have the correct aspect ratio. Figure [7b thus re presents a sl ice through 



the structure of the nebula in the plane along our line-of-sight (see ISmithl 120061 . for a more 
detailed explanation in the case of the bipolar lobes of Eta Carinae). This yields a tilt of 
the polar axis from the plane of the sky for Hen 3-401 of roughly 15°(±5°). In that case, the 
lateral expansion speed of the lobes (the radial expansion speed perpendicular to the polar 
axis) is V cxp ~15.5(±1) km s . In the images, the polar lobes have a width of roughly 2'.'5, 
which corresponds to 3.7xlO n (-Dfc pc ) km, where Dk pc is the heliocentric distance in kpc to 
Hen 3-401. Thus, the age of the nebula would be only about 380 xDk pc yr if the expansion 
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has been ballistic. Assuming a distance of 3 kpc (IBujarrabal fc Bachillerl Il99ll ) thus yields 
a nebular age of 1100 yr. 

There is some weak H 2 emission coming from the star as well. In fact, the very weak 
central peak appears to show some extended structure, with a possible tilt — i.e. it is 
slightly blueshifted toward the east and redshifted toward the west, as shown in Figure [8] in 
an expanded view of the region around the central star. This may be emission from an inner 
ring or disk surrounding the central star. This same type of tilted structure from the central 
H2 peak is also seen in similar spectra of the young PN M 2-9. In fact, it is remarkable 
that all aspects of the overall kinematic structure of H2 emission from Hen 3-401 in panel 
b of Figure [7J are almost identical to those of M 2-9 (jSmith et al. 2005 ) . Even the lateral 
expansion velocity of the lobes is roughly the same as for M 2-9. Although the apparent size 
of Hen 3-401 is smaller by a factor of ~4.5, it distance is larger by about the same amount, 
so the two objects have the same size, nebular age, and overall appearance. A detailed 
comparison of their central stars might therefore be quite illuminating. 

At most positions, H 2 emission from the far sides of the lobes is fainter than on the near 
side, especially where the far side is behind an area of bright emission on the near side. A 
simila r effect is seen in M2- 9, where it is probably due in part to extinction by dust in the 
lobes (jSmith fc Gehrzll2005l ). However, an asymmetry does appear in the brightness of the 
east-west lobes of Hen 3-401, with the brighter peak from the near side closer to the star in 
the west lobe, position « 2" while in the east lobe it is at position « 3"; this is in accord 
with what is seen in the H 2 image (Section 3.2). 

Slices in velocity across the positions +3" and —3" with width / .'26 are displayed in 
Figure Velocity widths of the arcs are narrower along the axis (Figure [7b ) than in the 
offset positions. This is due to projection effects in a Hubble flow. For a hollow cylinder 
with walls of constant thickness, the path length through the walls reaches a minimum for 
a sight line passing through the polar axis, translating to the smallest spread in velocity. 



(V 



These observations reveal a radial velocity for the system of Vh e i(H2) = —24 km s 1 

15 km s _1 . These are in good agreement with the 



LSR= 



=—36 km s 



with V, 



cxp 



millimeter-wave CO J=l— emiss ion-line val ues of Vlsr(CO) = — 3 6 km s 1 with V, 



16 km s 1 by lLoup et al.l (jl990l ) (although IBujarrabal fc Bachillerl (Il99ll ) measured the 



cxp 



somewhat different value of Vlsr(CO) 
.7=1-0 and CO .7=2-1 lines). 



28 km s 1 with V exp = 15 km s 1 using the CO 
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4-2.3. Rob 22 

Figure [TOl shows the long-slit Phoenix spectra of Rob 22, with one slit positioned along 
the middle of the lobes through the center of the nebula (panel b) and offset slit positions 
on either side. The very weak H 2 emission and the complex continuum structure in the 
long-slit spectra indicate that even a narrowband F212N image is actually dominated by 
dust structures seen in reflected continuum light. This is true over most of the nebula, with 
the exceptions being the N and S ends of the polar lobes. In the center slit position (b) there 
appear to be two continuum sources separated by ~0'.'6, with the northern one the brighter, 
while in the slit position to the west there is one continuum source. These are consistent 
with the H 2 image. These continuum sources likely represent light from the obscured central 
star scattered toward us by regions of enhanced density. The three panels show a slight 
positional offset in the emission pattern, in agreement with the "S" shape of the nebula. 

The spectra of the nebula show very faint H 2 emission from two hollow expanding lobes 
with an ellipsoidal velocity structure. The southern lobe is brighter and larger in H 2 , in agree- 
ment with the H 2 image. The nebula appears to have a polar axis very close to the plane of 
the sky, although the southern lobe is slightly blue-shifted and thus tilted slightly toward us. 
The lateral expansion speed of the lobes (difference in velocity between the front and back of 
the S lobe) is approximately 30 km s _1 . Thus V exp ~ 15 km s _1 and Vh e i(H 2 ) « +5 km s _1 , 
which transforms to Vlsr ~ — 7 km s _1 . These are in general agre ement with the OH maser 



measurement of Vlsr (OH) = — 6±3 km s 1 , V exp =20 km s 1 (lAUen. Hyland. &: Caswell 
1980l : ISilva et~aHll993h . However, they differ from t he molecular-line CO 7=1— measure- 



ment of V LS r(CO)=— 0.1 km s 1 , V exp ~ 35 km s 1 feujarrabal fc Bachiller 1991 ). although 



the investigators commented on the severe galactic CO contamination in this direction. 



5. DISCUSSION 

5.1. IRAS 16594-4656 

The PV diagrams for IRAS 16594—4656 indicate an expanding ellipsoidal nebula with 
its major axis nearly in the plane of the sky. This is a different orientation and morphology 
than had been concluded based upon the visual images of the nebula. In the visual images 
(see Fig. ITTk). the star appears bright and the SW lobe(s) appears brighter and perhaps 
larger than the NE lobe(s). This is in contrast to nebulae like AFGL 2688 (Egg nebula), 
IRAS 17150-3224, IRAS 17441-2411, and Rob 22, in which the star is obscured and the 
two lobes look more symmetrical in size, suggesting that the polar axis is very nearly in 
the plane of the sky. Inspection of the visual images had previously led us to conclude that 



-14- 



IRAS 16594—4656 was inclined at an intermediate orientation (~30°) wit h respect to th e 
plane of the sky with the SW lobe oriented in the direction of the observer (ISu et al.ll200ll ). 
Even the new H 2 image, with the SW lobe appearing larger than the NE lobe, suggests this, 
assuming symmetry in the lobes. 

However, other more recent observations at other wavelengths support our conclusion 
based upon kinematics that the axis of the main bipolar lobes, as delineated in the 



imag e, is nearly in the plane of the sky. A polarization study of this source (jUeta et al. 



20051 ) indicates a hollow shell with reflection from the walls. The polarization structure 
is aligned with the major axis of the nebula, with low polarization perpendicular to this, 
suggesting the presence of a dusty torus. A high-resolutio n mid-infrared study clearly shows 
the limb-brightened ends of such a torus (jVoIk et al.ll2006l ). The mid-infrared image does not 
show the elliptical structure expected of a torus seen at an intermediate angle, but rather 
shows a bar with bright ends, consistent with a torus seen nearly edge-on. The bright mid- 
IR emission is located at the two ends of the pinched waist seen in the H2 image, showing 
the presence of enhanced dust density which shapes the outflow. The approximately similar 
amounts of polarizations seen at opposite sides of the nebula are consistent with a nebula in 
which the bipolar axis is close to the plane of the sky rather than at an interme diate angle, 



and the authors suggest an angle of 15° with respect to the plane of the sky (jUeta et al. 



20051 ). Thus both the polarization and mid-infrared images are consistent with the result 
of this kinematic study, and the corresponding one of IVan de Steene et al.l (120081 ). that the 
bipolar lobes are aligned nearly in the plane of the sky. 

The H2 image shows limb-brightened lobes, consistent with the idea that they were 
carved out by collimated outflows. Based on the angular diameter across the lobes of 2'.'0 
and the measured V exp = 9 km s -1 , one can calculate an approximate expansion age for the 
lobes of 530 x/} fcpc yr. Using the polar velocity determined using SHAPE, along with the 
angular diameter between the two ends of the bipolar lobes of 4'.'0, leads to a expansion age 
of 740 xDkpc yr This leads to an estimated age for the lobes of 1400-1900 yr, using a distan ce 
of 2.6 kpc based on an assumed total luminosity of 10 4 L Q ( IHrivnak. Volk. fc Kwokl|2000l ). 



Comparing the morphology of the nebula in different wavelength regions shows several 
striking features. In visible light (Fig. [TTa). the nebula shows a basic bipolar structure, but 
with a series of three pairs of point-symmetric thin lobes. Thus at this wavelength it has 
been classified as point-symmetric and mu ltipolar. Portions of several circumstellar rings 
are also seen (IHrivnak. Kwok. fc Sull200ll ). In the J band (Fig. [TTb) the thin lobes are 
present and a portion of at least one ring is seen faintly. Seen more prominently, however, 
is a somewhat elliptical morphology extending from the northwest to the southeast that is 
also present in the visible image. However, it is not actually an ellipse, since the axes are 
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not separated by 180°; the northwest lobe is at a p osition angle o f 305° and the southeast 
one at 110°. The published H-band image is similar (ISu et al.ll2003l ). with the star relatively 
brighter compared to the fainter thin lobes. The H 2 image (Fig. [TTb). in contrast, shows a 
clear bipolar, "peanut" morphology with shorter, more spherical lobes. The major axis is at 
a position angle of 75° with the west lobe larger. Some of the structures that define these 
peanut-shaped lobes can be seen upon close examination in the V and H images. Among 
these is the bright "cap" at the end of the eastern lobe. As mentioned earlier, the mid-IR 
images at 10 and 20 /mi reveal two bright regions just outside the pinched waist of the 
"peanut." The mi d-IR emission also faintly delineates the edges of the two peanut-shaped 
lobes seen in H 2 . IVolk et al.l (120061 ). in comparing the mid-IR images to this H 2 image, 
concluded that the coincidence of the faint mid-IR emission with the H2 emission from the 
walls of these lobes suggested that the mid-IR emission was likely partially sh ock heated. 



The walls of these peanut-shaped lobes are clearly outlined in polarized light (jUeta et al. 
20071 ). agreeing with the position of the dust seen in the mid-IR. 



In a composite image (Fig. [TTH). we see the enhanced H2 emission from the walls 
of the peanut-shaped bipolar lobes and the H2 at the ends of the long, thin lobes. The 
non-uniformity in the emission from the walls may indicate lower density regions, through 
which jets may have carved holes from which the long, thin lobes may have emerged. This 
composite image also shows clearly several different structures seen in this object at different 
wavelengths, reinforcing the idea that the morphology of this nebula is more complex than 
simply symmetric bipolar lobes. It is more complex than either of the other PPNs of this 
study. 



In their excellent study of this object, IVan de Steene et al.l (120081 ) used their long-slit 
H 2 spectra together with our spectra (from the archive) to produce channel maps and from 
these to determine a more detailed three-dimensional image of the shell structure of IRAS 
16594—4656. While it basically consists of two hollow bipolar lobes joined at the waist, it is 
more complex, with some holes and extended structures (see their Fig. 8). These correspond 
to the brightness variations and the "cap" seen in our H 2 image. They also analyzed [OI] and 
[CI] emission lines, and making the reasonable assumption that these lines were produced in 
the same shock as the H 2 lines, they used these as diagnostics of the physical conditions in the 
shock region. From these, they deduce that the emission arises in very high density regions, 
3 x 10 6 < n R < 5 x 10 7 cm~ 3 . At these conditions of density, temperature (T ro t~14:00 K; 



Van de Steene &: van Hooill2003r). and velocity (8 km s -1 ), the excitation can be due to either 
C-type or J-type shocks (IVan de Steene et al.l 120081 ) . They also made long-slit observations 
of the [Fell] 1.645 /mi emission line, which indicate that this emission arises, not in the lobes, 
but from a wind close to the star. 
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5.2. Hen 3-401 



This nebula is very highly collimated, with the lobes extending ~14'.'5 from the central 
star with a length to width ratio for each lobe of ~7. The lobes in this object do not appear 
to be closed, but rather have "tattered" ends. The PV diagrams show a Hubble flow in the 
velocity in the lobes. This can arise du e to an eruptive ejecti on (ballistic motions) or (self- 
similar) hydrodynamic growth patterns (IBalick fc Frankll2002l ) . The velocities show that the 
eastern lobe is tilted away from us and the western lobe toward us. From the complementary 
image and spectroscopy, an estimate can be made of the velocity at the ends of the lobes 
and the inclination angle of the bipolar axis. With the measured l ength of the open- ended 
lobes of 14'.'5 and the width of 2.' 2 (from the sharpened V image, ISahai et al.lll999bl ). and 
the lateral lobe expansion velocity of 15.5 km s _1 , one can calculate a minimum expansion 
velocity at the ends of the lobes of ~200 km s _1 . The PV diagram along the bipolar axis 
(Fig. [7b) shows that the velocity at the ends of the two lobes differ by 20 km s _1 from 
that of the star. This implies an inclination angle of < 6° to the plane of the sky. This is 
somewhat smaller than we found from the appearance of the PV diagram in Section I4.2.2| 
but both clearly indicate a small inclination angle. 

In visible light, a "skirt" structure is seen surrounding the inner part s of the bipolar lobe s 
and extending coaxially with the lobes for ~3" from the star (discussed by lSahai et al.lll999bl ). 
This skirt also appears in the H-band image but not in H 2 . This is seen in Figure IT21 The 
wester n lobe is the brighte r in scattered light, with brightness ratios in the different bands of 
2.4 (V; ISahai et al.lll999bl ). 1.4 (J), and 1.4 (2.15 /mi). If one assumes an intrinsic symmetry 
in the two lobes, then this suggests some obscuration by an extended equatorial disk. In 
contrast, the eastern lobe is brighter in H 2 , with a ratio of 1.2. This may be a result of the 
differin g local shock condit ions in the nebula. There is a suggestion of a t orus from the optical 



1995 



image (ISahai et all Il999bt) and from the polarization measurements (jScarrott fc Scarrott 



Ueta et al.ll2007l ). A recent high-resolution mid-infrared imaging study shows that the 



mid-IR emission arises mostly from a 1" c ore, with weaker emission arising from the walls 
extending out to about 3" from the star (IMuthumariappan. Kwok. fc Volkl 120061 ) . These 
authors suggest that the core contains a flared disk, based upon the temperature gradient 
detected in the core. 

In the H band, the star is more prominent than in the visible and the inner region of 
the lobes is slightly less prominent. In visible light, there is a bright inner edge to the lobes, 
with enhanced light out to the ends of the skirt. This is also seen in the H band but not the 
H2 image. A faint halo is seen in both the visible and H-band images. The H 2 emission is 
particularly bright 2.0—5.5" from the star. These different regions of brightness are seen in 
the composite image (Fig. IT2ti). 
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Based upon the extremely large intensity ratio of the H2 1—0 S(l) to 2—1 S(l) line of 26 
( iGarcfa- Hernandez et al.1 120021 ). it is clear that the excitation mechanism is due to shocks. 
These presumably arise as the fast wind, that carved out the narrow lobes, produces shocks 
i n the walls that excite th e H 2 . Ha wings are observed to extend out to ~1600 km s _1 
( jGarcia-Lario et al.lll999bl ). 

The young PN M 2-9 is very simi l ar to H en 3-401 in visible appearance, as has been 
remarked upon by iGarcfa-Lario et al.1 (Il999bl ), and also in its mid -infrared appearance 
(although with not quite as pinched a waist; ISmith fc Gehra 120051 ) and its kinematics. 
It possesses smooth, almost cylindrica l lobes, a bright central star, and a thin H2 zone 
( IHora fc Latterlll994i ; ISmith et al.ll2005l ). A recent study of winds driven by magnetic pres- 
sure from a toroidal field was able to reproduce well the morphology of these two objects 
( IGarcia-Segura et al.ll2005l ). 



5.3. Rob 22 



The PV diagrams for Rob 22 show an expanding nebula with the H2 emission arising 
from the ends of the lobes and the walls. The bipolar axis lies nearly in the plane of the 
sky with the south lobe expanding slightly toward us compared with the north lobe. Given 
the expansion velocity of ~15 km s" 1 and the angular size of the lobes of ~2.0", this leads 
to a kinematic age of ~320 xDk pc yr. For a distance of ~2 kpc (jAUen. Hyland. fc Caswell 
19801 ). this yields an age for the lobes of ~640 yr. 



In the visual image (Figure [T3"a ; ISahai et al.lll999al ). the south lobe is brighter than 
the north lobe with a ratio of 1.6. In the near-infrared, the north lobe is slightly brighter, 
with ratios of 1.05 in H (F160W; Figure [13b), 1.1 in 2.15 /im, and 1.2 in the H2 image 
(Figure [T3b). If the reduced relative brightness in the south lobe in the V band is due to 
extinction by the dense equatorial disk, then the disk would need to be extremely large if it 
were orthogonal to the axis of the outflow, which as we have seen is very close to the plane 
of the sky. It may be more likely that there is simply an intrinsic asymmetry in the lobes 
that contributes most of the difference. 

If we compare the different images (Fig. [T3|) . we s ee a strong point symmetry in the "S" 
shape. Even at 7.7 and 8.7 fim this "S" shape is seen (IMuthu. Kwok. fc Volkll2004j ) . Such a 
pattern can be formed by collimati on of an outflow by a precessing disk. The Ha emission 
has wings to at least ±450 km s -1 (jAUen. Hyland. fc Caswelllll980l ). 



The V image shows a distinct dark band separating the lobes. In the H image, this dark 
lane is not so prominent, and a bright, compact clump appears. This bright clump is red 
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but not point-like, and likely represents warm dust heated by the obscured central star. The 
nebula is brighter close to the star, with a particularly bright region NE of the center of the 
nebula. A similar morphology to the H-band image i s seen in a recent ad aptive optics very 
high resolution (R~0.06") K-band image obtained by lLagadec et al.l (120081 ); this shows some 
finer detail and in particular a bright linear extension toward the SE from the very bright 
region near the center, defining the southern part of the "V" -shaped central obscured region. 
The halo is evi dent in both the V and H images. A faint region in the N lobe (referred to 
as a "spur" by ISahai et al.l (jl999al )) is seen in both the V and H images, although in the V 
image it is outlined by emission farther out in the lobe while in the H image it is outlined 
by emission closer to the center. In the H 2 image, the emission is particularly strong in the 
same region NE of the center, with a less bright region on the other side (SW) of the center. 
In addition, the H 2 is particularly bright along the ends of the lobes. This can arise due 
to shock excitation produced as the collimated wind, that may have produced the S-shaped 
nebula, interacts with the AGB wind of the halo. Polarization measurements sho w scattering 



from the ends of the lobes, in addition to the presence of a equatorial dust torus (jUeta et al. 
2007h . 



The faint region "spur" in the N lobe is evident in all three images. It is interesting 
that in the OH observations (IDyer et al.ll200ll ). the more negative velocities arise from this 
faint region; perhaps it represents an evacuated hole punched through the nebula. 



T he OH emission appears to be confine d primarily to the equatorial regions (ISahai et al. 



1999al : lDver et al]l2001 



Zijlstra et al.ll200ll ). where it may represent the expansion or rotation 



of a disk. The mixed C and O circumstellar chemistry might also suggest a disk, as this 
explanation has been advocated as the cau se of the mixed ch emistry in several other post- 
AGB objects, such as the Red Rectangle (Ivan Winckell 120031 ) . The detection of hot dust, 
manifested in the spectral energy distribution of this object (but not in the other two), is 
consistent with a disk. Or it may be that this object is in transition between an O-rich 
chemistry to a C-rich one, with the PAH emission arising closer to the star, as may be the 
case for PN with [WC] central stars. 



5.4. Comparison With Other Spatial-Kinematic Studies of PPNs 

A small number of other bipolar PPNs have also been studied spatially or kinematically 
in H 2 . We discuss these for comparison. 



AFGL 2688 ("Egg" nebula) has been imaged using HST-NICMOS (ISahai et al.lll998f ) 
and shows H 2 emission from clumps in its spindle-shaped closed lobes, especially near their 



-19- 



outer ends. This is consistent with the emission arising along the walls of the lobes. It 
also shows strong e mission from an extended equatorial region. A kinematic study by 
Kastner et all (120011 ) shows that the highest velocity H 2 along the polar axis resides close 
to the star and that the velocity of the H 2 decreases along the polar axis away from the 
star. They find a similar effect for AFGL 618. They interpret this as resulting from a faster 
collimated wind colliding with slower-moving AG B material. T h e equ atorial H2 shows one 
side to be blue- and the other red-shifted, which iKastner et all (120011 ) interpret indicating 
a combination of expansion a nd rotation. An alternate interpretation would be multiple 
bipolar jets along the equator (ICox et al.ll2003l ) rather than rotation, and this seems to be 
the more likely explanation. This decrease with distance along the polar axis seen in AFGL 
2688 and AFGL 618 is in contrast with what is seen in Hen 3-401. Hen 3-401 displays a clear 
trend of the H 2 velocity increasing in velocity with distance away from the star, in accord 
with what one would expect from an unhindered ejection of material from the region of the 
central star. 



Davis et all (120051 ) recently made spatially-resolved high-resolution (R~20,000 or 16 



km s _1 ) H 2 observations of four bipolar PPNs. They observed them with the spectrograph 
slit oriented along the bipolar axis and then offset slightly to each side, similar to our 
observing pattern for Hen 3—401 and Rob 22. In two of these objects, M 1—2 and IRAS 
17441—2411, one sees a distinct velocity difference between the two lobes (AV of ~55 and 
~10 km s^ 1 , respectively), for IRAS 17243—1755 H 2 is detected only weakly from one lobe, 
and for IRAS 17150—3224 the two lobes have similar velocities (bipolar axis in the plane of 
the sky). For this last object, they also suggested that there may be rotation in the lobes 
based on systematic velocity differences seen in the offset positions. 



Hrivnak et al.1 (120061 ) also studied IRAS 17150—3224 with higher spectral resolution, 



and they obtained a high-resolution (0'.'2) i/ST-NICMOS H 2 image. The image shows the 
H 2 emission to come primarily from clumps near the ends of the lobes, but these are not 
resolved in the spectra, which basically show the integrated emission and velocity of the 
clumps. The two lobes have a similar velocity, as in the above study, with slight velocity 
differences in the offset positions due to small effects of the different H 2 clumps rather than 
rotation. 



As we have referred to previously, IVan de Steene et al.1 (120081 ) also recently studied the 
kinematics of IRAS 16594—4656. They also observed with the Phoenix spectrograph on 
Gemini-South at even somewhat higher spectral and spatial resolution than we did. Their 
spectra show a similar picture of an expanding H 2 ellipsoid with V exp = 8 km s _1 . They 
combined their spectra with ours to produce a detailed three-dimensional image of the shell 
structure. 
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In none of these earlier studies does one see such clearly-resolved spatial-kinematic struc- 
tures as are seen in the new Gemini-South Phoenix data of IRAS 16594—4656, Hen 3—401, 
and Rob 22. In all three the emission from the walls of the lobes is clearly traced in the PV 
diagrams. The small slit size and high spectral resolution, together with the high spatial 
resolution images, result in exceptional kinematic detail for a PPN. 



6. SUMMARY AND CONCLUSIONS 

The high-resolution HST-NICMOS H2 images of these three PPNs have allowed us to de- 
termine the spatial location of their H 2 emission and the spatially-resolved H 2 spectroscopic 
observations have allowed us to determine the kinematics of the H 2 emitting regions. In all 
three cases, the systemic and expansion velocities are similar to those determined from the 
molecular-line CO or OH measurements. We find the following results for these three PPNs. 

IRAS 16594—4-656: While the V and H-band images show a complex multi-lobe struc- 
ture, the H 2 images shows a clear bipolar, peanut shape , although with variations in density 



holes") and structure (see IVan de Steene et al.ll2008l ). The H 2 emission originates along 



the walls of the lobes (sides and ends), with fainter emission from more distant (ejected?) 
clumps. The PV diagram shows the H 2 to arise in an expanding ellipsoidal velocity struc- 
ture, which is in contrast to the bilobes of the density structure. The kinematics indicates 
that the bipolar lobes are nearly in the plane of the sky (if»10°); this differs from the earlier 
interpretations of the lobes as being at some intermediate orientation, but it is consistent 
with recent mid-IR imaging and near-IR polarization studies. The lobes are estimated to 
have an age of ~1600 yr. 

Hen 3-401: The H 2 emission originates from the sides of lobes, which have open ends. 
The lobes are tilted somewhat to the plane of the sky (~10— 15°), with the western lobe 
moving toward us, and they show an increasing velocity with radial distance (Hubble flow). 
An estimated age for the lobes is ~1100 yr. The open ends and unhindered outflow may be 
the consequence of its higher degree of collimation (greater linear momentum). 

Rob 22: The H 2 emission originate primarily from ends of the S shaped nebula and from 
regions of the S shape near the obscured central star. The nebula is nearly in the plane of 
the sky, consistent with the absence of a visible star due to obscuration by a disk. An age 
of ~630 yr is estimated for the lobes. 

H 2 surveys of PPNs have shown that H 2 emission is commonly found in those with a 
bipolar morphology. As can be seen from this study, the combination of H 2 high-resolution 
images and spatially-resolved, high-resolution spectra provides valuable insight into the 
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structure and shaping mechanisms for these bipolar nebulae. 
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Table 1. HST NICMOS Observing Log 



IRAS ID Other Names Observation Total Exposure Times (s) a 

Date F212N F215N F110W F160W F222M 



10178-5958 Hen 3-401 1998 Mar 6 1728 1728 ■ • • 576 1152 
10197-5750 Rob 22 1998 May 18 1344 1344 • • • 320 768 
16594-4656 • • • 2002 Aug 13 1215 1215 192 



''Note that in tw o of our previous papers involving NICMOS images ( ISu et al.l 12003 ; 



Hrivnak et al.ll2006l ). we incorrectly listed as the total exposure times values that were one- 



third of the correct values. Thus for our F212N and F215N images of IRAS 17150-3224 , for 
example, the total exposure times were 2448 s, not 816 s as listed (IHrivnak et al.ll2006l ). 
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Table 2. Photometric Calibration Values a 



Filter 


A c 


AA 


PHOTFLAM 


PHOTFNU 


ZP(Vega) 


PSF FWHM b 




fxm 


fim 


erg cm~ 2 / um _1 DN _1 


Jy sec DN" 1 


Jy 


a 


F212N-new 


2.1213 


0.0087 


1.9135E-14 


2.8722E-5 


664.7 


0.19 


F215N-new 


2.1487 


0.0079 


2.1031E-14 


3.2389E-5 


645.1 


0.20 


F212N-old 


2.1213 


0.0088 


2.4272E-14 


3.6432E-5 


664.7 


0.20 


F215N-old 


2.1487 


0.0079 


2.5533E-14 


3.9323E-5 


645.1 


0.20 


FllOW-new 


1.1235 


0.1630 


2.9059E-15 


1.2234E-6 


1784.4 


0.12 


F160W-old 


1.6060 


0.1177 


2.3956E-15 


2.0610E-6 


1040.7 


0.14 



a The photometric calibration values were taken from the online NICMOS photometry perfor- 
mance for post-NICMOS (new) cooling System installation (operating temperature of detector 
~77K) and pre-NICMOS (old) Cooling System installation (operating temperature of detector 
~62K). 

b Based upon a Gaussian fit to the observed field stars. 
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Table 3. High- Resolution H 2 Spectroscopic Observing Log 



Object 




Date 


P.A. (°) a 


Position 


^exp (s) 


No. Obs. 


Rob 22 




2003 Jan 14 


17 


through star 


500 


2 


Rob 22 




2003 Jan 14 


17 


0"7 W of star 


500 


2 


Rob 22 




2003 Jan 14 


17 


0f!7 E of star 


500 


2 


Hen 3-401 




2003 Dec 17 


73 


through star 


500 


2 


Hen 3-401 




2003 Dec 17 


73 


Qf!7 S of star 


500 


2 


Hen 3-401 




2003 Dec 17 


73 


0"7 N of star 


500 


2 


Hen 3-401 




2003 Dec 17 


163 


through star 


500 


2 


IRAS 16594- 


-4656 


2003 May 08 


52 (b) 


through star 


600 


2 


IRAS 16594- 


-4656 


2003 May 08 


33 (c) 


through star 


600 


2 


IRAS 16594- 


-4656 


2003 May 08 


345 (e) 


0"9 W of star 


600 


2 


IRAS 16594- 


-4656 


2003 May 15 


345 (d) 


through star 


600 


3 


IRAS 16594- 


-4656 


2003 May 15 


72 (a) 


through star 


600 


2 



a The letters in parenthesis for IRAS 16594—4656 refer to the slit orientations as displayed 
in Figure HI 

b Exposure times per observation. 
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Fig. 1— The flux-calibrated F212N (left) and F215N (right) images of IRAS 16594-4656 
(top), Hen 3-401 (middle), and Rob 22 (bottom). The images are all in reversed, logarithmic 
scale. 
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Wavelength (n-m) 
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1 > 1 1 i 1 ' 1 ' 1 1 1 1 < 1 1 1 1 1 1 > 1 1 1 1 1 > i 1 1 1 1 

2.1 2.12 2.14 2.16 2.18 

Wavelength (\un) 



Fig. 2. — The medium- resolutio n spectra of IRAS 16594—4656 , Rob 22, and Hen 3-401 
( Van de Steene &: van Hooi 2003 ; Garcfa-Hernandez et al.l [2OO2I ) . showing the H 2 and the 
Br7 emission features, along with the HST NICMOS F212N and F215N filter profiles. The 
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Fig. 3. — In the left-hand panels are shown the resulting flux-calibrated H 2 images of IRAS 
16594—4656 (top), Hen 3-401 (center), and Rob 22 (bottom). In the right-hand panel of 
each is shown the H2 image with the position of the Phoenix slits overlaid. The labelling of 
the slits corresponds to the labelling in the position- velocity diagrams (Figs. 141 [71 [TUj) . 
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Fig. 4. — The H 2 position-velocity diagram of IRAS 16594—4656. Panels a - e correspond 
to the five slit positions shown in Fig. [31 
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Fig. 5. — PV diagram for IRAS 16594—4656. On the left side is shown the observed data, 
in the middle is shown the model results with a constant velocity plus slight Hubble flow ( V 
= (8 + 2.4r/l") km s _1 ), and on the right is shown the results with a pure Hubble flow ( V 
= 12r/l" km s _1 ). The top row is for slit position a and the bottom row for slit postion d 
in Fig. [3j Note the good agreement of the shape of the the middle panels (constant velocity 
with slight Hubble flow) with the observed data. (To fit the brightness would involve a 
modification in the surface emissivity from the constant value assumed.) The model results 
have been convolved with the seeing and velocity resolution to produce the suitable images 
to compare with the observations. 
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Fig. 6. — Superimposed slices in velocity at position along slit = 0" across the H 2 position- 
velocity diagram of IRAS 16594—4656. (The differences in flux levels likely arise from dif- 
ferences in seeing and transparency.) 
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Fig. 7. — The H2 position- velocity diagram of Hen 3-401. Panels a - c correspond to the 
three slit positions shown in Fig. [3j 
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Fig. 8. — Enlarged H2 position- velocity diagram of the region around the central star of Hen 
3-401. 
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Fig. 9. — Slices in velocity at slit positions +3" and —3" from the central star and width 
0.26" across the H 2 position- velocity diagram of Hen 3-401; see Fig. [7] for reference. 
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Fig. 10. — The H 2 position-velocity diagram of Rob 22. Panels a - c correspond to the three 
slit positions shown in Fig. [31 Note that the gray-scale in each panel has been adjusted 
arbitrarily based upon the strongest emission levels; since the eastern offset position (a) does 
not contain a strong continuum source, the lower level emission in the lobes is enhanced in 
that panel relative to the others. 
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Fig. 11. — Images of IRAS 16594—4656 at various wavelengths: (a) F606W (~Fband; A e // 
= 0.603 jum , < AA > = 0.150 fim ), (b) F110W (~#band), and (c) H 2 emission (from Fig. 
3). Figure (d) is a color composite of the three, with blue representing the Fband, green 
the J band, and red the H2 emission. To better show the details in the central region, it is 
shown as an inset with the intensity scaled down by a factor of ten and then displayed with 
a logarithmic scale. 
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Fig. 12. — Images of Hen 3—401 at various wavelengths: (a) F606W (~ Fband; A e // = 0.603 
fim , < AA > = 0.150 /im ), (b) F160W (~i/band), and (c) H 2 emission (from Fig. 3). 
Figure (d) is a color composite of the three, with blue representing the Fband, green the H 
band, and red the H 2 emission. 
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Fig. 13. — Images of Rob 22 at various wavelengths: (a) F606W (~V A band; A e jj = 0.603 
/iin , < AA > = 0.150 fim ), (b) F160W band), and (c) H 2 emission (from Fig. 3). 
Figure (d) is a color composite of the three, with blue representing the Vband, green the H 
band, and red the H 2 emission. 



